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ABSTRACT 

We report the results of experiments aimed at reducing the major problem with cooling flow models of 
rich cluster X-ray sources: the fact that most of the cooled gas or its products have not been found. Here 
we show that much of the X-ray emission usually attributed to cooling flows can, in fact, be modeled by 
a power-law component which is indicative of a source(s) other than thermal bremsstrahlung from the 
intracluster medium. We find that adequate simultaneous fits to ROSAT PSPCB and ASCA GIS/SIS 
spectra of the central regions of ten clusters are obtained for two-component models that includes a 
thermal plasma component that is attributable to hot intracluster gas and a power-law component that 
is likely generated by compact sources and/or extended non-thermal emission. For five of the clusters that 
purportedly have massive cooling flows, the best-fit models have power-law components that contribute 
~ 30 % of the total flux (0.14 - 10.0 keV) within the central 3 arcminutes. Because cooling flow mass 
deposition rates are inferred from X-ray fluxes, our finding opens the possibility of significantly reducing 
cooling rates. 

Subject headings: cooling flows — X-rays: galaxies — X-rays: galaxies: clusters 



1. INTRODUCTION 

It is well established that the diffuse X-ray emission from 
clusters of galaxies implies the existence of a hot (10 7-8 
K), thermally radiating intracluster medium (ICM). In the 
dense centers of many clusters this gas is expected to cool 
radiatively in less than a Hubble time. As pressure sup- 
port decreases the gas will gravitate towards the bottom of 
the cluster potential well, a phenomenon that has come to 
be known as a "cooling flow" (cf. Fabian 1994). Such cool- 
ing flows can be detected by direct X-ray imaging, where 
they are seen as a surface brightness enhancement in the 
X-ray emission from a region typically ~ 100 - 200 kpc in 
radius, and by spectroscopic observations which reveal gas 
at cooler temperatures than the surrounding ICM. The ex- 
cess X-ray luminosity in the cooling flow region is typically 
10 42 < L coo i < 10 45 ergs s _1 , corresponding to roughly 10 
- 50 % of the total cluster X-ray emission (Peres et al. 
1998). The mass deposition rate is M oc L coo ;T _1 , where 
T is the gas temperature at r coo ;, the cooling radius. Esti- 
mated cooling flow rates in the centers of rich clusters vary 
from 30 - 300 M yr" 1 (White, Jones, & Forman 1997) 
and in some cases reach mass deposition rates of 1000 M Q 
yr _1 or more (Allen 2000). 

A long-standing problem with the cooling flow scenario 
has been the fact that the large quantities of cool gas ex- 
pected to accumulate from the flow have rarely been de- 
tected. Searches for cool gas in atomic or molecular form 
have generally proven unsuccessful (e.g., Heckman et al. 
1989; Jaffe 1990; McNamara, Bregman, & O'Connell 1990; 
O'Dea et al. 1994; Voit & Donahue 1995; Jaffe & Bremer 
1997; Donahue et al. 2000), with reported upper limits 
two orders of magnitude smaller than the masses implied 
by X-ray mass deposition rates (although see Edge 2001). 
While some support for the cooling flow picture is pro- 



vided by low energy X-ray spectra of clusters which ap- 
pear to require larger absorption than can be accounted 
for by foreground gas alone (White et al. 1991; Allen et al. 
1993; Irwin & Sarazin 1995), it now appears that the frac- 
tion of clusters that actually possess such excess absorp- 
tion may be much less than previously thought (Arabad- 
jis & Bregman 2000, hereafter AB2000; Bohringer et al. 
2001a). Furthermore, recent high- resolution X-ray spectra 
obtained with the XMM-Newton and Chandra satellites 
have failed to show expected emission lines from cooled gas 
at temperatures below a few keV (Fabian et al. 2001; Pe- 
terson et al. 2001; Tamura et al. 2001) and it has also been 
found that local isothermality matches the spectral data of 
"cooling flow" regions better than the continually declin- 
ing temperature profile predicted by the standard cooling 
flow model (e.g. Molendi & Pizzolato 2001; Bohringer et 
al. 2001a). A variety of mechanisms have been proposed 
to account for the lack of detectable cool gas in cooling 
flows; these include reheating the gas, mixing, differen- 
tial absorption, efficient conversion of cooling flow gas into 
low mass stars, inhomogeneous metallicity distributions, 
or disruption of cooling flows by recent subcluster mergers 
(see recent discussions in Fabian et al. 2001 and Peterson 
et al. 2001). As of yet, it is unclear whether any of these 
mechanisms are capable of resolving the aforementioned 
problems with the cooling flow model. 

One characteristic that is common to the mechanisms 
mentioned above is that they have mainly been introduced 
in an "after the fact" manner. That is, they assume the ob- 
served mass deposition rates are basically correct but that 
something is happening to the gas after (or while) it cools, 
be it reheating or conversion into low mass stars or some 
other process. Perhaps a simpler way of addressing the 
problems of the cooling flow model are to look for ways to 
reduce the mass deposition rates from the beginning. One 
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way to achieve this is to reduce the excess emission {Lcooi^) 
associated with cooling flows by attributing some (or all) 
of it to X-ray sources other than thermal bremsstrahlung 
or line emission from the ICM. But are there sources lumi- 
nous enough to reproduce L coo i , can they be found within 
galaxy clusters, and, most importantly, can they mimic 
hot cluster gas? 

Soon after X-ray emission was first detected from galaxy 
clusters, Katz (1976) proposed that much of this emis- 
sion might come from individual compact sources such 
as those known to exist in our own Galaxy. However, 
the subsequent discovery of iron emission lines in X-ray 
spectra (Forman & Jones 1982, and references therein) 
clearly demonstrated that thermal bremsstrahlung from 
a diffuse intracluster gas is the dominant source of X- 
ray emission in clusters, and interest in compact sources 
quickly waned. Nevertheless, compact sources must surely 
contribute some fraction of the total X-ray emission from 
galaxy clusters. A growing body of evidence indicates that 
a substantial intergalactic population of stars, star clusters 
and stellar remnants is present in the cores of rich clusters, 
and some of these objects will be X-ray sources. For ex- 
ample, Ferguson, Tanvir, & von Hippel (1998) showed the 
existence of intergalactic stars in the Virgo cluster, and 
planetary nebulae associated with this population have 
also been detected (e.g., Arnaboldi et al. 1996; Ciardullo 
et al. 1998). Similarly, West et al. (1995) hypothesized 
the existence of a large population of intergalactic glob- 
ular clusters in the cores of galaxy clusters based on ob- 
served galaxy-to-galaxy variations in globular cluster pop- 
ulations. Furthermore, there is mounting evidence that 
many cluster galaxies have been destroyed over a Hub- 
ble time, victims of collisions with other galaxies or dis- 
ruption by the mean tidal field of the cluster potential 
in which they reside (e.g., Thompson & Gregory 1993; 
Lopez-Cruz et al. 1997; Gregg & West 1998; Moore et al. 
1999; Calcaneo- Rodin et al. 2000). As galaxies are dis- 
rupted they will spill their contents throughout intraclus- 
ter space, including X-ray sources such as X-ray binaries, 
globular clusters, and black holes. Such disruption is likely 
to be most frequent in the dense cluster core, and hence 
would be expected to produce a substantial pool of intra- 
cluster material, including compact X-ray sources, at the 
cluster center. Any significant contribution by these com- 
pact sources to the X-ray flux will have to be taken into 
account when estimating the properties of ICM, including 
temperatures, densities, and cooling flow mass deposition 
rates. 

An intracluster population of compact objects would 
have a different spectral signature than thermal emission 
from intracluster gas. For instance, the spectra of large 
spiral and elliptical galaxies reveal the presence of hard 
X-ray emission which is thought to be produced by either 
high mass X-ray binaries or accretion onto massive black 
holes and is adequately represented by power-laws having 
T - 0.6 - 3.0 (where F v cx v~ T ; Colbert & Mushotzky 
1999; Allen, Di Matteo, & Fabian 2000). Luminous AGN, 
on the other hand, are often fit by power-laws with Y ~ 
1.6 - 2.0 (Nandra et al. 1997). Indeed, the spectral signa- 
tures of such objects have already been identified in several 
X-ray clusters by Allen et al. (2001). Using ASCA data, 
these authors found very luminous sources with power-law 



spectral indices characteristic of compact sources. 

The recent discovery of luminous extended non-thermal 
X-ray emission in several weak (or non-) cooling flow clus- 
ters (Rephaeli, Gruber, & Blanco 1999; Fusco-Fcmiano et 
al. 1999; 2000; 2001) suggests that one should also look for 
such emission in massive cooling flow clusters. At present, 
the origin of the diffuse non-thermal emission is still a 
matter of some debate. Suggested mechanisms include 
inverse-Compton scattering of CMB photons up to X-ray 
energies via a population of highly relativistic electrons, 
and non-thermal bremsstrahlung, likewise, by a popula- 
tion of relativistic electrons (Sarazin 2001; Rephaeli 2001). 
Those highly relativistic electrons may have been accel- 
erated by recent subcluster merger events which shocked 
the ICM or may have been ejected from a central massive 
black hole. The discovery of X-ray 'holes' at the posi- 
tions of radio lobes of a central compact object in several 
massive cooling flow clusters (e.g. Bohringcr et al. 1993, 
McNamara et al. 2000; Fabian et al. 2000) certainly lends 
credence to the second hypothesis. Recently, Harris et 
al. (2000) (Chandra data) and Bohringer et al. (2001b) 
(XMM-Newton data) have even identified extended non- 
thermal emission directly associated with radio jets/lobes 
in 3C295 and M87, respectively. 

The spectral signature of extended non-thermal emis- 
sion is also different from that of thermal emission. To 
date, the extended non-thermal components of X-ray spec- 
tra of clusters have adequately been modeled by power- 
laws having r ~ 0.5 — 6.0 (Fusco-Femiano et al. 1999; 
2000; 2001). As with luminous compact sources, a proper 
deduction of cluster properties should explicitly include 
the effects of extended non-thermal emission in the anal- 
ysis. 

In the present work we assess the case for a power-law 
component to the X-ray spectra of cluster centers, in par- 
ticular those with cooling flows. We show that simulta- 
neous fits of ROSAT PSPCB and ASCA GIS/SIS spec- 
tra, which have been used to infer the presence of cooling 
flows in many clusters to date, are in fact consistent with a 
significant contribution from compact sources and/or ex- 
tended non-thermal emission. Using ROSAT HRI data, we 
have placed limits on the contribution of any bright point 
sources within the central regions of these clusters and 
conclude that the power-law components are most likely 
produced by either a collection of many low-luminosity 
compact sources (with Lx < 5 x 10 42 ergs s _1 ) or by dif- 
fuse non-thermal emission or by some combination of the 
two. We shall assume H a = 50 km s _1 Mpc -1 and q = 
throughout. 

2. DATA SELECTION AND REDUCTION 

Ten well-studied X-ray luminous galaxy clusters (Table 
1) were chosen to span a wide range of cooling flow mass 
deposition rates. Five of these have essentially no cooling 
flow (NCF), while the other five cooling flow (CF) clusters 
have calculated mass deposition rates ranging from 189 to 
645 solar masses per year. Observations by the ROSAT 
Position Sensitive Proportional Counter (PSPCB) and 
High Resolution Imager (HRI) and the ASCA Gas Imaging 
Spectrometer (GIS) and Solid-state Imaging Spectrometer 
(SIS) were obtained from the HEASARC public archives 
and reduced using FTOOLS v5.01. The combination of 
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ROSAT and ASCA data provides good sensitivity to both 
soft and hard X-ray emission, spanning the total energy 
range 0.14 - 10.0 keV, and should provide a strong test of 
power-law components in the X-ray spectra of clusters. 

2.1. ROSAT observations 

The ROSAT PSPCB data were corrected for spatial and 
temporal gain fluctuations in the detector. X-ray spec- 
tra, spanning the energy range 0.14 - 2.04 keV, of all ten 
clusters were extracted from circular regions centered on 
the peak of the emission with a constant extraction radius 
(R) of 3 arcminutes. Background spectra were selected 
from source-free circular regions of radius R. All spectra 
were corrected for the effects of cosmic rays in the detec- 
tor (Snowden et al. 1992), vignetting, spatial variations 
in detector efficiency, and deadtime. The background- 
subtracted spectra were then grouped to ensure a mini- 
mum of 20 counts per radial bin so that reliable x 2 statis- 
tics could be computed. 

The ROSAT HRI data was reduced in a manner simi- 
lar to the PSPCB data. The HRI data was used to place 
limits on the contribution of any point sources within the 
central regions of these clusters. Unfortunately, the HRI 
lacks spectral resolution and was not used in the main 
spectral analysis. A description of the analysis of HRI 
data is given in Section 5. No HRI data were available for 
A2063 or A3532. 

2.2. ASCA observations 

The screened events from "revision 2" processing of 
ASCA CIS/SIS data were used. Additional cleaning pro- 
cedures (discussed in the ASCA ABC Guide) were imple- 
mented, including gain and deadtime corrections. For the 
CIS detectors (GIS2 and GIS3), X-ray spectra spanning 
the energy range 0.7 - 10.0 keV were extracted from the 
same region of sky as that of the ROSAT PSPCB source 
spectra. For the SIS detectors (SIS0 and SIS1), X-ray 
spectra spanning the range 2.0 - 10.0 keV were extracted, 
thus avoiding complications with the degrading efficiency 
of the SIS detectors at low energies. Background spectra 
for both the SIS and CIS detectors were extracted from 
the "blank sky" observations made during the performance 
testing stage of the mission. Spectra were extracted from 
the same regions of the detectors as that of the source 
spectra (thus minimizing position-dependent effects asso- 
ciated with background subtraction). The background- 
subtracted spectra were also grouped (or binned) to en- 
sure x 2 statistics could be used. For the clusters A2142 
and A3667 SIS1 data were ignored because of an excess of 
hot and flickering pixels. 

3. SPECTRAL ANALYSIS 

Reduced ROSAT PSPCB and ASCA GIS/SIS spectra 
were fit simultaneously using models within the software 
package XSPEC vll (Arnaud 1996). Each cluster was 
fit with a variety of models based on MEKAL thermal 
plasma routine (Mewe & Gronenschild 1985; Kaastra 1992; 
Liedahl et al. 1995), and/or redshifted power-laws. In ad- 
dition, we used the popular cooling flow model comprising 
a thermal plasma (MEKAL) and cooling flow component 
implemented by the MKCFLOW routine within XSPEC 
(Mushotzky & Szymkowiak 1988). Fits were obtained by 



varying the temperature (kT) and metallicity (Z) of ther- 
mal plasma component (s) and by varying the photon index 
(r) of any power-law components. Photoelectric absorp- 
tion (Balucinska-Church & McCammon 1992) was applied 
to the spectral models through two separate methods: i) 
by applying only one absorption component to the entire 
spectral model and fixing the HI column density (Nh) at 
the known Galactic value (see Table 1) and then letting it 
vary to see if the fit could be improved; ii) by applying one 
absorption component to the entire spectral model (with 
Nh fixed at the Galactic value) and a separate redshifted 
variable absorption component (Nh : pl) to any power-law 
components via the ZPHABS model in XSPEC. The sec- 
ond method allows for intrinsic absorption to be associated 
with any power-law components. Spectra of nearby high 
mass X-ray binaries and accreting massive black holes pro- 
vide strong evidence for heavy intrinsic absorption, with 
Nh ~ 10 21 ~ 23 cm~ 2 above the Galactic column (Colbert 
& Mushotzky 1999; Iwasawa et al. 2000; Akylas et al. 
2001). We therefore regard the addition of a free intrin- 
sic absorption parameter as an important feature of our 
spectral analysis. 

Henceforth we will use TP, PL, and CF to refer, re- 
spectively, to thermal plasma (MEKAL), power-law, and 
MKCFLOW components. For example, in this nota- 
tion TPPL identifies a model having thermal plasma and 
power-law components. Goodness of fit was judged by 
the reduced chi-square (xV) and by visual inspection of 
residual plots. When additional components were added 
to a spectral model, for example adding a power-law com- 
ponent to a thermal plasma component to make TPPL, 
the corresponding F-statistic was calculated in order to 
judge whether the addition was reasonable. Unless stated 
otherwise, the derived relative contributions of each of the 
components to the total flux refer to the energy range 0.14 
- 10.0 keV stated, whereas the luminosities presented re- 
fer to the ASCA GIS band only (0.7 - 10.0 keV), since 
the spectra were fit simultaneously and not combined. All 
luminosities presented were corrected for Galactic absorp- 
tion only. 

To check our data reduction procedures, spectra from 
each of the detectors were also fit individually and no sig- 
nificant deviations were found. We have also attempted 
to duplicate published results by Arabadjis & Bregman 
(2000) for three clusters derived using the same data 
(ROSAT PSPCB) and the same software. The results, 
which refer to the inner 3 arcminute radius, are presented 
in Table 2. Here, the hotter TP component is fixed at the 
White, Jones, & Forman (1997) value, Z = 0.5Z Q , N H 
is fixed at the Stark et al. (1992) value, and the energy 
range is 0.14 - 2.4 keV and 0.14 - 2.04 keV for AB2000 and 
our study, respectively. The comparison shows excellent 
agreement, with derived temperatures and mass deposi- 
tion rates having the same values to within one standard 
deviation. 



4. RESULTS 

4.1. Non-Cooling Flow Clusters 

We find that single-temperature TP models provide ac- 
ceptable fits to the spectra of A119, A2063, A3266, A3532, 
and A3667 (1.00 < xl < 1-17) with temperatures in very 
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good agreement with previous studies (White, Jones, & 
Forman 1997; Markevitch et al. 1998). Significant im- 
provement is not usually obtained by varying the hydrogen 
column density from the known Galactic value reported 
by Dickey & Lockman (1990) or Stark ct al. (1992). The 
derived metallicities range from 0.2 < Z < 0.4 which is 
typical of the intracluster medium (Allen & Fabian 1997). 

In general, both pure PL models with or without vari- 
able intrinsic absorption do not provide acceptable fits to 
the NCF cluster spectra with Xr si 1-4- This result is not 
unexpected since it leaves no room for thermal emission 
from the ICM. With intrinsic absorption included, an ac- 
ceptable fit is obtained for A3532 with x 2 ~ 1-05, however. 
This cluster is located at the center of the Shapley Con- 
centration and is thought to be undergoing a merger with 
A3528. The diffuse non-thermal emission generated from 
such a merger might possibly account for the ability of the 
simple PL model to fit the spectrum of A3532. Unfor- 
tunately, no HRI data exists for this cluster so we cannot 
place limits on the contribution of any point sources within 
the central 3 arcminutes. 
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Fig.l. The best fitting TPPL model on the spectrum of A3667. The 
power-law component has a photon index of 1.30 and contributes ks 
22 % of the flux over the range 0.14 - 10.0 kcV. Only the ROSAT 
PSPCB data are presented here. 

In the best fitting TPPL models, the PL component typ- 
ically contributes 1 - 14% (mean of 7 %) when absorbed by 
the Galactic column only or 10 - 25 % (mean of 20 %) when 
intrinsic absorption is included (excluding A3532 which 
has an anomalously large PL component). We note that 
in virtually every case the inclusion of intrinsic absorption 
noticeably improves the quality of the fit, making it supe- 
rior to even the simple TP models (see Table 3). A mean 
F-statistic of « 10 (> 95 % confidence) was found when the 
PL component with intrinsic absorption was included in 
the fit indicating that its addition is, indeed, reasonable. 
The corresponding luminosities of intrinsically absorbed 
PL components range from Lx, 0.7-10.0 ~ 5 x 10 42 — 7x 10 43 
erg s" 1 with a mean value of L x ,o.7-w.o ~ 2.2 x 10 43 
erg s -1 (corrected for Galactic absorption only). Pho- 
ton indices for these components are typical of compact 
accretion-driven X-ray sources and extended non-thermal 
emission. Interestingly enough, the amount intrinsic ab- 
sorption ranges from 1.8 x 10 21 < N h ,pl < 1-4 x 10 22 
cm~ 2 which is typical of optically thick compact sources. 
We also note that the temperature of the ICM tends to 



decrease significantly when a power-law component is in- 
cluded. This trend was previously reported by Allen et 
al. (2001) when including a power-law component in their 
spectral models. Figure 1 presents the best fit TPPL 
model to the spectrum of A3667. For clarity, only the 
ROSAT PSPCB data are plotted. 

As expected, standard cooling flow (TPCF) models of 
our non-cooling flow sample yield best fits with M = M Q 
yr~ 4 to within one sigma. The mass deposition rate for 
the modest cooling flow cluster A2063, however, is 25 ± 
12 Mq yr" 1 . The cooler component in our TPTP models 
typically contributes zero flux to within one sigma (~ 10 % 
for A2063 & A3667). Neither of these models provided fits 
significantly better than those of the standard TP models. 

In summary, we confirm earlier investigations which find 
that TP models provide satisfactory fits to ROSAT and 
ASCA spectra of the centers of non-cooling flow clusters. 
However, TPPL models which include intrinsic absorption 
generally provide fits of superior quality, and predict that 
the power-law component contributes 20 percent of the 
total flux. The photon indices of this component are typi- 
cal of accreting X-ray sources and/or diffuse non-thermal 
X-ray emission which usually range from 1.0 < T < 3.0. 

4.2. Cooling Flow Clusters 

We have fit the spectra of the CF clusters with the same 
variety of models used in the previous section, but also 
with an additional model; TPTPPL. This model repre- 
sents the emission from the observed region as the sum of 
the hot ICM, a diminished cooling flow (second TP com- 
ponent), and intermixed compact sources and/or extended 
non-thermal emission characterized by the PL component. 

Simple TP models generally provide fairly reasonable 
fits to the spectra with 1.1 < x 2 < 1.25. Derived tem- 
peratures and metallicities are in good agreement with 
previous studies (White, Jones, & Forman 1997, Allen & 
Fabian 1997; Sarazin & McNamara 1997; Markevitch ct 
al. 1998; Allen 2000; AB2000, Allen et al. 2001). None 
of the fits were significantly improved by letting the value 
of Nh vary. Hence, we find no evidence for excess absorp- 
tion in these clusters, which is also consistent with several 
other investigations (Sarazin 1996; Allen & Fabian 1997; 
Sarazin & McNamara 1997; Allen 2000; AB2000). 

As anticipated, pure PL models (with or without intrin- 
sic absorption) do not provide acceptable fits to the X-ray 
spectra of any of the CF clusters. Typically, Xr si 1-5 f° r 
this model. 

Our TPPL models fit the spectra of all 5 CF clusters 
very well, with 1.0 < xl < 1-2 (see Table 4). With the ex- 
ception of A2029, which is best fit by simple TP models, 
the PL component typically contributes 5 - 30 % (mean of 
20 %) of the total flux when absorbed only by the Galactic 
column. The inclusion of intrinsic absorption significantly 
improves the quality of the fit and raises the contribution 
of the PL component to 9 - 47 % (mean of 27 %). The 
mean F-statistic value for the addition of this power-law 
component is « 50 which suggests it is significant (> 99 % 
confidence). The luminosity (0.7 - 10.0 kcV) of PL compo- 
nents in these four clusters span w 0.9 — 8 x 10 44 ergs s _1 , 
with a mean value of « 3.8 x 10 44 ergs s _1 . The photon 
indices of PL components are reasonably well constrained 
and range from 1.0 < T < 2.4 (mean of 1.8) (see Table 
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4). The implied intrinsic column densities acting on the 
PL component are rather high for three of the CF clusters 
(A644, A1689, & A2597) and are consistent with high- 
mass compact sources. As with the NCF clusters, a lower 
temperature is found for the thermal component than is 
found from the simple TP models. 

TPTP models were also implemented. The quality of 
the fits provided by these models is generally comparable 
to that of TPPL models. Calculation of the F-statistic 
suggests that the addition of the extra TP component is 
significant. Thus, we confirm previous reports that multi- 
phase thermal plasmas fit the spectra of CF clusters much 
better than that of isothermal models. The temperatures 
and metallicites of the two thermal components are consis- 
tent with the results of previous studies (White, Jones, & 
Forman 1997, Markevitch et al. 1998; Allen 2000; AB2000; 
Allen ct al. 2001). For the cooler component in particu- 
lar, we find 1.4 < kT < 2.2 and a mean contribution of 
roughly 22 % of the total flux, which is typical of cooling 
flow clusters (Peres et al. 1998). 

The standard cooling flow model (TPCF) provides ad- 
equate fits to the spectra of the CF clusters, as expected. 
Derived mass deposition rates, temperatures, and metallic- 
ities are generally consistent with previous studies (White, 
Jones, & Forman 1997; Peres et al. 1998; Allen et al. 
1999; Allen 2000; AB2000). The best fitting TPCF model 
to the spectrum of A644 is presented in Figure 2. Only 
the ASCA SIS0 and GIS2 data are plotted. A compari- 
son of the TPCF and TPPL models is presented in Table 
4. It is clear that the TPPL models fit the spectra of the 
5 CF clusters just as well as standard TPCF models and 
without any additional complexity 
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Fig. 2. ASCA SIS0 and GIS2 spectra of the CF cluster A644. The 
SIS0 spectrum is from 2.0 - 10.0 keV while the GIS2 spectrum ranges 
from 0.7 - 10.0 keV. Both data sets are shown in the residuals in the 
bottom part of the plot. The best fitting TPCF model predicts a 
mass deposition rate of 200 Mq yr~L 

Finally, the TPTPPL model also fits the X-ray spectra 
of these clusters well, with 1.0 < Xr ^ 1-2- When absorbed 
by the Galactic column only, the PL component typically 
contributes 2 - 5 % of the total flux inside 3 arcminutes. 
The addition of intrinsic absorption on the PL component, 
again, improves the fit and the mean contribution of the 
PL component is raised to w 32 %. When compared to 
the TPPL model, the mean F-statistic value is about 5 (« 
68 % confidence) which suggests the addition of a second 



TP component only improves the fit marginally. The lu- 
minosity (0.7 - 10.0 keV) of the PL component ranges from 
1 - 6 x 10 44 ergs s _1 with a mean value of w 2.2 x 10 44 
ergs s _1 . The power-law components have photon indices 
of 1.8 < r < 3.3 with a mean of 2.3. With the exception 
of A2029 (which has virtually no contribution from the PL 
component), the intrinsic column densities on the PL com- 
ponent are > 10 21 cm~ 2 . The temperature of the cooler 
TP component did not change significantly from the value 
derived from TPTP models. However, it's contribution to 
the flux was generally lower than that of TPTP models 
(mean - 10 %). 

In summary, we find that in all cases the TPPL and 
TPTPPL models fit the ROSAT and ASCA spectra of the 
five CF clusters just as well as the standard cooling flow 
model (TPCF) or the two-temperature thermal plasma 
(TPTP). The PL components of these models are quite 
luminous (contributing between 27 - 32 % of the total 
flux) and characterized by photon indices typical of com- 
pact accreting X-ray sources and/or extended non-thermal 
emission. This, then, is the main conclusion of our work: 
ROSAT and ASCA spectra can be satisfactorily modeled by 
attributing a significant part of the luminosity within the 
central regions of CF clusters to sources other than cooling 
gas. 

5. DISCUSSION 

Our analysis shows that luminous (Lx. 0.7-10.0 ~ 10 42-44 
ergs s _1 ) power-law components with 1.0 < T < 3.0 might 
be present in the centers of rich clusters. Including intrin- 
sic absorption dramatically improves the quality of the fit; 
in these circumstances, there is an apparent dichotomy 
in the importance of the PL components between NCF 
and CF clusters. Not only is the typical PL component 
in the CF cluster more luminous by about a factor of 10 
than that of its NCF counterpart, but it also contributes 
a larger fraction of the total X-ray flux within the central 
3 arcminutes of its cluster (between 27 - 32 % for CF clus- 
ters as opposed to 20 % for NCF clusters). This may be 
consistent with the fact that the cores of CF clusters are 
more often observed to be active (i.e. radio loud), perhaps 
implying the presence of a central massive black hole(s) 
and a surrounding diffuse plasma of relativistic electrons, 
than the cores of NCF clusters (Burns 1990). It is also 
worth noting that the improvement of the fit to the spec- 
tra of CF clusters far exceeds the improvement of the fit to 
the spectra of NCF clusters when a power-law component 
is included in the spectral model. 

The most attractive feature of the TPPL models is that 
one does not expect large quantities of atomic or molecular 
gas to be present in the centers of clusters, nor does one 
expect to see emission lines below a few keV, because cool- 
ing is not significant. This is consistent with current ob- 
servations. On the other hand, evidence which apparently 
supports the existence of cooling flows (e.g. excess absorp- 
tion, detection of CO line emission in extreme CF clusters) 
may contradict this model. However, the success of the 
single-temperature plasma + power-law model in fitting 
the spectra of the central regions of rich clusters, regardless 
of cooling flow status, has also been reported by a number 
of authors using a variety of X-ray instruments, including 
Markevitch & Vikhlinin (1997) (ASCA and ROSAT data), 
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Markevitch et al. (1998) (ASCA data), Rephacli, Gru- 
ber & Blanco (1999) (RXTE data), Guainazzi & Molendi 
(1999) (Beppo-SAX data), Fusco-Fcmiano ct al. (1999; 
2000; 2001) (Beppo-SAX data), and recently by Sambruna 
et al. (2000) who used Chandra to uncover a previously 
hidden AGN in Hydra A with Lx,2-io > 10 42 ergs s^ 1 and 
T <~ 1.7. However, our study is the first that systemati- 
cally addresses the problems of the standard cooling flow 
model by implementing the single-temperature plasma + 
power-law model. 

The TPTPPL models, which might represent a com- 
promise between the standard cooling flow model and our 
favoured TPPL models, also provide good fits to the spec- 
tra of the CF clusters. However, the addition of the sec- 
ond, cooler TP component only marginally improves the fit 
and at the expense of adding much more complexity (and 
uncertainty) to the model. Furthermore, the predicted 
temperatures and relative contributions of the cooler TP 
component in the TPTPPL models suggest that emission 
lines below a few keV should be present, however they have 
failed to turn up in Chandra and XMM spectra. Finally, 
because the mean contribution of the cooler component is 
- 10 % , as compared with 20 % in the TPTP models, 
this implies that mass deposition rates are only roughly 
reduced by a factor of two. Using models similar to our 
TPTPPL models, Allen et al. (2001) reached the same 
conclusion for a number of clusters in their sample. These 
reduced deposition rates are still far too large to be con- 
sistent with the small amounts of atomic and molecular 
gas found in cooling flow clusters. We expect that Chan- 
dra and XMM-Newton will place much tighter constraints 
on the parameters of the complex TPTPPL models and it 
will be interesting to see if the models can be reconciled 
with problems of the standard cooling flow model. 

Assuming that a power-law component is, indeed, 
present in the spectra of these clusters, some of the possi- 
ble candidates are: 

• Seyfert-like AGN could contribute most, if not all, 
of the PL component luminosity found in our study. 
These objects can be quite luminous and are of- 
ten observed to have Lx ~ 10 41 ~ 44 ergs s _1 with 
power-law indices of 1.6 - 2.0 (Nandra et al. 1997). 
In fact, Allen et al. 2001 have suggested that AGN 
are the likely candidate for the observed PL com- 
ponent in their sample. Likewise, the luminosities 
and steepness of the PL components in our sam- 
ple are generally consistent with emission from a 
central AGN(s). 

• Diffuse non-thermal emission, which has possibly 
been seen in the central regions of Coma, A2256, 
and A2199. Searches of those clusters using Beppo- 
SAX and RXTE have yielded very luminous (w 
10 42-44 ergs s _1 ) sources which are described by 
power-laws with T ~ 0.5 — 6 and may possibly be 
linked to power-law components observed in the 
ASCA band (Rephaeli, Gruber, & Blanco 1999; 
Fusco-Fcmiano et al. 1999; 2000; 2001). Unfortu- 
nately, X-ray spectra usually have to be extracted 
from regions much larger than that of the central 
galaxy and, thus, it is difficult to draw conclu- 
sions about the source of the power-law component. 



A VLA survey of more than 200 nearby clusters 
found that only 29 had extended radio emission on 
Mpc scales (Giovannini, Tordi, & Feretti 1999) per- 
haps implying that extended non-thermal emission 
is not that common in galaxy clusters. It should be 
noted, however, that the detection of true intraclus- 
ter radio emission is exceptionally difficult, as one 
must remove bright radio galaxies from the analy- 
sis (Rephaeli 2001). Furthermore, the non-thermal 
emission may be associated with the radio galaxies 
themselves (e.g. in the radio lobes). 

• Massive black holes associated with elliptical and 
lenticular galaxies. The regions sampled in our 
study generally include a single giant elliptical or 
cD galaxy plus a number of large galaxies, proba- 
bly E or SO galaxies. Colbert & Mushotzky (1999) 
and Allen, Di Matteo, & Fabian (2000) have re- 
cently discovered hard power-law components in 
several nearby gE galaxies with 10 37 < Lx < 10 42 
ergs s _1 with photon indices of 0.5 - 3.0. Such 
sources are typically less luminous than Seyfert 
nuclei, perhaps as a result of low radiative effi- 
ciency from an advection-dominated accretion flow 
(ADAF; Narayan & Yi). While several of these 
sources could, in principle, account for the PL com- 
ponent in the NCF clusters they would have to be 
scaled up by a factor of 100 to account for the emis- 
sion in the CF clusters. Faint galaxies, on the other 
hand, are much more numerous and many of them 
might harbor less massive central black holes. The 
Local Group elliptical M32, whose central object 
has Lx ~ 10 37 erg s" 1 (Loewenstein et al. 1998) of- 
fers a possible guide to the X-ray luminosity of faint 
cluster members. Hence ~ 10 4 — 10 6 such galaxies, 
each with a supermassivc black hole, would be re- 
quired to account for the power-law component of 
the clusters studied here. On this basis, faint galax- 
ies seem incapable of producing significant power- 
law emission inside the regions we have observed. 
We note, however, that if massive black holes orig- 
inally belonging to bright or faint galaxies are now 
scattered throughout intracluster space (i.e. no 
longer associated with optically visible galaxies) as 
the result of past mergers or stripping both of these 
limits may be relaxed. 

• Galactic X-ray binaries, having T « — 2.5 (White 
et al. 1988, Verbunt et al. 1995), and which in our 
galaxy are most commonly associated with globular 
clusters. Intergalactic globular clusters, if similar 
to globular clusters in the Milky Way, would have 
to be present in large numbers to collectively ac- 
count for ~ 10 42 ~ 44 erg s _1 attributable to the PL 
component in our models. Assuming a typical in- 
tergalactic globular contains a binary X-ray source 
with L x ~ 10 35 ergs s" 1 (Hertz & Wood 1985) and 
the specific frequency of globulars in the central gE 
galaxy is Sx = 10 (Harris 1991) then even a very 
bright central galaxy with = -24 would have 
a total X-ray luminosity of <~ 10 39 ergs s -1 , too 
low by a factor of 10 4-5 . Even the most optimistic 
models of intracluster globular clusters of West et 
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al. (1995) predict ~ 10 4 — 10 5 such clusters, signifi- 
cantly less than is required to account for the entire 
PL component. 

Hence, it appears that there are probably two or three 
sources capable of producing the luminosity of the PL com- 
ponents of our models: massive black holes associated (or 
once associated) with large elliptical galaxies, diffuse non- 
thermal emission, and AGN. Does existing data favor one 
type of source over the others? 

To determine the luminosity of any point sources within 
the central 3 arcminutes we used ROSAT HRI data. 
Two different procedures were implemented: I) Spectra 
of clearly visible point sources were extracted from a 5" 
region centered on the source and grouped into one large 
energy bin (0.5 - 2.0 keV). The spectra were then fit within 
the XSPEC package using a simple power-law model with 
T = 1.8 and absorbed by the Galactic column only. The 
0.5 - 2.0 keV flux was then converted into a 0.7 - 10.0 keV 
unabsorbed flux using the program W3PIMMS. If more 
than one point source was present the individual fluxes 
were summed. As an example, Figure 3 is the HRI ob- 
servation of A3667, smoothed with a 5" Gaussian beam. 
There are at least two bright compact sources within the 
central 3 arcminutes of this cluster.; II) To place limits 
on the contribution of point sources which are hidden in 
the bright diffuse background (what we refer to as "barely 
visible" point sources) we calculated the mean standard 
deviation of counts per pixel (within 3') and then scaled 
this to a standard deviation per circular region of diam- 
eter = 5" (comparable to the HRI psf). The latter value 
was expressed as unabsorbed flux (0.7 - 10.0 keV) using 
the program W3PIMMS and assuming a simple power-law 
spectral model with T = 1.8 absorbed only by the Galac- 
tic column listed in Table 1. We regard the result as a 
1-sigma limit to the contribution of any point sources. 




Fig. 3. ROSAT HRI observation of A3667 (centered on R.A. 
20 h 12 m 31 s , Dec. -56°49T2", J2000). Two extremely bright com- 
pact sources are visible within the central 3' but together are roughly 
a factor of 7 less luminous then the spectrally determined power-law 
component in the best fitting TPPL model. 



The results of HRI data analysis are presented in Ta- 



ble 5. One or more point sources are easily visible within 
the central 3 arcminutes of all three NCF clusters having 
HRI observations. It is clear from a comparison with the 
luminosities of the observed PL component from the best 
fitting TPPL models (with intrinsic absorption), however, 
that the bright point sources cannot be solely responsible 
for the power-law emission, as they contribute only « 5 % 
of the PL component luminosity. Significant contributions 
are required from low- luminosity point sources and/or dif- 
fuse non-thermal emission. The third column in Table 5 
(Method II) presents the 3-sigma limit of "barely visible" 
point sources, suggesting that at least 10 2 ~ 3 of them are 
needed (in addition to the visible point sources) to account 
for the entire power-law component. 

There are no visible point sources within the central 3' of 
the 5 CF clusters studied. This is almost certainly the re- 
sult of the CF clusters being much more centrally peaked 
than the NCF clusters. At least 10 2-3 "barely visible" 
point sources are required to account for the observed PL 
component in these clusters. 

Perhaps another way to differentiate between compact 
sources and diffuse non-thermal emission is through the 
presence or absence of flux variability, since most luminous 
accreting X-ray sources are known to vary significantly on 
short time scales (e.g. Nandra et al. 1997). However, be- 
cause the PL component typically contributes only « 20 
% this might make detection of variability difficult in clus- 
ters (especially in the CF clusters). Furthermore, if the PL 
component is produced by 10 2 ~ 3 point sources, and the 
variations in each are random, it will make the variation 
in the collective luminosity quite small, perhaps rendering 
this method unusable. 

If compact or diffuse non-thermal sources provide signif- 
icant emission from the centers of rich clusters, it will be 
necessary to revise a number of conclusions derived from 
X-ray data. Fully exploring the consequences is beyond 
the scope of the present paper; however, we will mention a 
few important issues. Allen (2000) reports that a deprojec- 
tion analysis of ROSAT HRI data yields cooling flow mass 
deposition rates that agree reasonably well with those de- 
rived from ASCA spectra. A consistent re-analysis would 
require first subtracting the power-law contribution from 
the X-ray images, but of course the spatial distribution 
of the power-law emission must first be understood. The 
presence of emission lines from heavy ions originally led 
astronomers to accept thermal emission as the dominant 
mechanism in rich cluster ICMs. Accretion-driven X-ray 
sources do exhibit line emission (Kallman, Vrtilek, & Kahn 
1989; Narayan & Raymond 1999), but the lines are pro- 
duced under quite different physical conditions than those 
within the canonical hot ICM. The implications for ICM 
abundance estimates and enrichment scenarios will require 
careful study. Likewise, previously derived temperatures 
of the ICM will also be affected (as demonstrated in Ta- 
ble 3). Might emission from compact sources account for 
the excess hard X-radiation recently discovered (Xu et al. 
1998; Ikebe et al. 1997) in the centers of several clusters? 
The next generation of X-ray detectors and satellites (i.e. 
XMM and Chandra) will be needed to map the X-ray emis- 
sion with adequate spatial resolution to address this and 
other questions raised by the present work. 
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6. CONCLUSIONS 

We have shown that ROSAT PSPCB and ASCA 
GIS2 / GIS3 spectra of ten X-ray luminous clusters permit 
a significant contribution from compact sources and/or dif- 
fuse non-thermal emission. Although this conclusion holds 
irrespective of the presence of a cooling flow, we find that 
power-law components are especially luminous in clusters 
with massive cooling flows. Our study indicates that in- 
ferred cooling flow mass deposition rates may be lowered 
substantially, thus ameliorating, or possibly resolving, the 
long-standing problem of the missing cool gas. Exist- 
ing X-ray observations of large elliptical galaxies, AGN, 
and extended non-thermal emission suggest that any of 
these sources could account for the observed luminosity of 
the power-law component. However, the absence of ex- 
tremely bright point sources in the ROSAT HRI observa- 
tions rules out the possibility that the power-law emission 



comes from a few extremely bright compact sources. More 
likely is that most of the power-law emission is produced by 
some combination of diffuse and low-luminosity compact 
sources, such as massive black holes. The deep imaging ca- 
pabilities and superior spectral resolution of Chandra and 
XMM should be able to place much tighter constraints on 
the relative contributions of power-law sources and ther- 
mal bremsstrahlung to cluster X-ray emission and, per- 
haps, even determine the source of the power-law emission. 
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Table 1 

Clusters investigated in this study. 



Cluster 


a 

Z 




Type 


M 


kl 


Notes 


A 1 1 Q 


0449 


^ 1 7 

O.LI 


IN V_.' JT 


n 
u 


o.o 




A644 


0.0704 


8.00 e 


CF 


189 


6.5 


minor merger event 9 


A1689 


0.1832 


1.81 


CF 


645 


8.7 




A2029 


0.0773 


3.14 


CF 


556 


7.4 


non-thermal emission not detected' 1 


A2063 


0.0353 


2.90 e 


NCF 


37 


4.1 


subcluster merger? 1 


A2142 


0.0909 


4.16 


CF 


350 


11.4 


merger? & hard X-ray emission detected* 


A2597 


0.0852 


2.49 


CF 


271 


9.1 


powerful cD radio galaxy' 


A3266 


0.0589 


3.00 e 


NCF 





6.2 


major merger" 1 


A3532 


0.0554 


5.98 


NCF 





4.7 


merger™ 


A3667 


0.0556 


4.00 e 


NCF 





7.1 


hard X-ray emission found 



a Rcdshifts from the NASA/IPAC Extragalactic database 

b Galactic HI column density (in units of 10 20 cm~ 2 ) from Dickey & Lockman (1990), except where noted. 
c Mass deposition rate (in M Q yr^ 1 ) from Peres et al. (1998) 

d Single-phase ICM temperature (in keV) reported by White, Jones, & Forman (1997). 

c Galactic HI column density (in units of 10 20 cm~ 2 ) from Stark ct al. (1992). 

'Feretti et al. (1999) 

SBaucr & Sarazin (2000) 

h Molendi & Dc Grandi (1999) 

'Oegerle & Hill (1994) 

JMarkevitch et al. (2000) 

k Bazzano et al. (1984) 

'Sarazin et al. (1995) 

m Flores, Quintana, & Way (2000) 

"Bardelli, Zucca, & Baldi (2001) 

°Fusco-Femiano ct al. (2001) 



Table 2 
Comparison with AB2000. 



Cluster 


Model 


kT coo i 


M b 


Xr 


kTcool 


M d 


Xr 


A119 


TP 






0.987 






1.011 




TPCF 


0.08 ± oo 


± 3 


0.994 


0.01 ± 151 


0.6 ± 5 


1.026 


A2029 


TPTP 


1.29 ± 0.15 




1.592 


1.19 ± 0.31 




1.499 




TPCF 


0.08 ± 6.97 


132 ± 31 


1.911 


0.03 ± 4.82 


170 ± 43 


1.772 


A2142 


TPCF 


0.08 ± 7.41 


188 ± 43 


1.190 


0.10 ± 9.16 


230 ± 73 


1.130 



a Temperature (in keV) of the cooler TP component found by AB2000. 
b Mass deposition rate (in M Q yr -1 ) found by AB2000. 
c Temperature (in keV) of the cooler TP component found in this study. 
d Mass deposition rate (in M Q yr -1 ) found in this study. 
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Table 3 

Comparison of the TP and TPPL models for the NCF clusters. 



Cluster 


Model 


kT a 


r b 




PL % c 


N h ,pl d 


9/7 c 

X /d-o.f. 


A119 


TP 


5.91 ±0.20 










467.4/446 




TPPL 


3.96 ±0.78 


1.55 ±0 


17 


25 % 


29.5 ± 17.6 


454.0/443 


A2063 


TP 


3.09 ±0.05 










945.6/807 




TPPL 


2.49 ±0.23 


1.41 ±0 


27 


8.5 % 


47.5 ±71.5 


933.3/804 


A3266 


TP 


8.48 ±0.24 










1160.5/1117 




TPPL 


4.23 ±0.86 


1.33 ±0 


12 


23 % 


140.0 ±35.2 


1158.0/1114 


A3532 


TP 


4.55 ±0.13 










776.4/746 




TPPL 


3.71 ± 1.14 


2.12 ±0 


10 


87.5 % 


18.0 ±2.70 


756.0/743 


A3667 


TP 


5.64 ±0.14 










809.0/791 




TPPL 


3.22 ±0.32 


1.30 ±0 


15 


22 % 


88.0 ±40.4 


797.6/788 



a Temperature (in keV) of the thermal plasma component. 
b Photon index of the power-law component. 

c Contribution of the PL component to the total flux (0.14 - 10.0 keV). 
d Intrinsic HI column density (in units of 10 20 cm~ 2 ) acting on PL component. 



Table 4 

Comparison of the TPPL and TPCF models for the CF clusters. 



Cluster 


Model 


M a 


r b 


PL % c 


N h ,pl d 


X 2 /d.o.f. 


A644 


TPPL 




2.37 ±0.78 


35 % 


24.1 ±20.6 


1855.7/1613 




TPCF 


200 ± 71 








1841.2/1613 


A1689 


TPPL 




1.02 ±0.29 


9 % 


166.0 ±71.5 


1643.6/1409 




TPCF 


504 ± 149 








1653.9/1409 


A2029 


TPPL 




2.13 ± 1.44 


5 % 


3.34 ±6.40 


1460.3/1149 




TPCF 


427 ± 60 








1462.1/1149 


A2142 


TPPL 




1.67 ±0.09 


47 % 


2.14 ±0.54 


1010.4/868 




TPCF 


307 ± 90 








1013.2/868 


A2597 


TPPL 




1.89 ±0.45 


17% 


177.0 ±72.0 


1305.1/1281 




TPCF 


217 ±23 








1302.3/1281 



a Mass deposition rate (in Mq yr 1 ). 
b Photon index of the power-law component. 

c Contribution of the PL component to the total flux (0.14 - 10.0 keV). 
d Intrinsic HI column density (in units of 10 20 cm~ 2 ) acting on PL component. 



Alternative to Cooling Flows? 



Table 5 

Limits in cgs units on point source contribution from HRI data. 



Cluster Method I a Method II b 



L 



X,PL 



A119 

A644 
A1689 
A2029 
A2142 
A2597 
A3266 
A3667 



8.12 x 10 41 



1.47 x 10 42 
8.01 x 10 42 



< 5.40 x 10 40 

< 2.99 x 10 41 

< 2.09 x 10 42 

< 5.34 x 10 41 

< 6.45 x 10 41 

< 7.02 x 10 41 

< 8.40 x 10 40 

< 8.97 x 10 40 



1.67 



+ 1.46 
-0.97 
~.74 



X 10 



43 



2.59I 2 ^ 4 x 10 44 
2.95t 4: ™ x 10 44 
8.65lf 6 f x 10 43 
8.8lt 4 ' 22 x 10 44 
O.MtUl x 10 44 



43 



6.28Z5 ii x 10' 

5-70±|i| x 10 43 



a Luminosity of visible point sources. 

b Uppcr limit of luminosity of "barely visible" point sources within 
the HRI psf. 

c Luminosity of the PL component of the best fitting TPPL models 
with intrinsic absorption. 



